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Probability distributions which emerge from the formalism of nonextensive statistical mechanics 
have been applied to a variety of problems. In this paper we unite modeling of such distributions with 
the model of widespread 1// noise. We propose a class of nonlinear stochastic differential equations 
giving both the g-exponential or g-Gaussian distributions of signal intensity, revealing long-range 
correlations and 1/ behavior of the power spectral density. The superstatistical framework to get 
1//^ noise with g-exponential and g-Gaussian distributions of the signal intensity in is proposed, as 
well. 
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I. INTRODUCTION 



Stationary stochastic processes and signals are prevalent across many fields of science and engineering. Many 
complex systems show large fluctuations of macroscopic quantities that follow non-Gaussian, heavy-tailed, power-law 
distributions with the power-law temporal correlations, scaling, and the (multi)fractal features The power- 

law distributions, scaling, self-similarity and fractality are sometimes related both with the nonextensive statistical 
mechanics and with the power-law behavior of the power spectral density, i.e., l/f^ noise (see, e.g., [sl l9l-[Tl|. 
and references herein). 

There exist a number of systems, involving long-range interactions, long-range memory and anomalous diffusion, 
that possess anomalous properties in view of traditional Boltzmann-Gibbs statistical mechanics. Nonextensive statis- 
tical mechanics represents a consistent theoretical background for the investigation of some properties, like fractality, 
multifractality, self-similarity, long-range dependencies and so on, of such complex systems, {6^-S\- Concepts related 
with nonextensive statistical mechanics have found applications in a variety of disciplines including physics, chemistry, 
biology, mathematics, economics, informatics, and the interdisciplinary field of complex systems (see, e.g., (12, - 14.] and 
references herein). 

The nonextensive statistical mechanics framework is based on the entropic form 

. _ i-gb(^)N^ 

where p{z) is the probability density of finding the system with the parameter z. Entropy ([T]) is an extension of the 
Boltzmann-Gibbs entropy Sbg — — p{z)^i^p{z)dz, which restores from Eq. ^ at q — I [6,?]. 

By applying the standard variational principle on entropy ^ with the constraints Jl^^ p{z)dz = 1 and 

z-[p{z)]^dz ^ 
J^-[p{z)]'^dz 

where is the generalized second-order moment [T5l - [l7j . one obtains the g-Gaussian distribution density 

p(z) = Aexp,(-Bz2). (3) 
Here expq(-) is the g-exponential function defined as 

e^Pqix) = [1 + {1 - q)x]]^ , (4) 

with [(. . .)]_!_ = (...) if (...)> 0, and zero otherwise. Asymptotically, as x — > oo, exp^(x) ~ x~^, where X = {q — 1)~^, 
i.e., we have the power-law distribution. The (more) generalized entropies and distribution functions are introduced 
in Refs. [11 [11. 

Statistics associated to Eqs. ([T|)-(j3]) has been successfully applied to phenomena with the scale-invariant geometry, 
like in low-dimensional dissipative and conservative maps (20l . [2l| . anomalous diffusion [2^ . turbulent flows [23| . 
Langevin dynamics with fluctuating temperature [23 - l26j . long-range many-body classical Hamiltonians 27 1, and to 
the financial systems [1^, [l^ . 

For the modeling of distributions of the nonextensive statistical mechanics, the nonlinear Fokker-Planck equations 
and corresponding nonlinear stochastic differential equations (SDEs) '55, 's^, SDEs with additive and multiplicative 
noises [Slj, [s^l, with multiplicative noise only [H, and with fluctuating friction forces [l^ have been proposed. 
However, the exhibition of the long-range correlations and 1//^ noise has not been observed. 

The phrases "1// noise", "1// fluctuations", and "flicker noise" refer to the phenomenon, having the power spectral 
density at low frequencies / of signals of the form S{f) ~ 1//*^, with (3 being a system-dependent parameter. Signals 
with 0.5 < P < 1.5 are found widely in nature, occurring in physics, electronics, astrophysics, geophysics, economics, 
biology, psychology, language and even music (33 - [38| (see also references in paper [l^). The case of /? = 1, or "pink 
noise", is the one of the most interesting. The widespread occurrence of processes exhibiting 1// noise suggests that 
a generic, at least mathematical explanation of such phenomena might exist. 

One common way for describing stochastic evolution and properties of complex systems, is by means of generalized 
stochastic differential equations of motion (39l - l4l| . These nondeterministic equations of motion are used in many 
systems of interest, such as simulating the Brownian motion in statistical mechanics, in fundamental aspects of 
synergetics and biological systems, field theory models, the financial systems, and in other areas [3. l39l. |42|. |43| . 



3 



The purpose of this paper is to model together both, the Tsalhs distributions and 1// noise, using the same 
nonUnear stochastic differential equations. The superstatistical approach for modeling of such processes is proposed, 
as well. 

We considered a class of nonlinear stochastic differential equations giving the power-law behavior of the probability 
density function (PDF) of the signal intensity and of the power spectral density (l//*^ noise) in any desirably wide 
range of frequency. Modifications these equations by introduction of an additional parameter yields Brownian-like 
motion for small values of the signal and avoids power-law divergence of the signal distribution, while preserving 1/ 
behavior of the power spectral density. The PDF of the signal generated by modified SDEs is q-exponential or q- 
Gaussian distribution of the nonextensive statistical mechanics. The superstatistical framework using a fast dynamics 
with the slowly changing parameter described by nonlinear stochastic differential equations can retain 1/ behavior 
of the power spectral density as well. When the PDF of the rapidly changing variable is exponential or Gaussian, we 
obtain g-exponential or g-Gaussian long-term stationary PDF of the signal, respectively. 



II. NONLINEAR STOCHASTIC DIFFERENTIAL EQUATION GENERATING ASYMPTOTICALLY 

POWER-LAW SIGNALS WITH 1//" NOISE 



Starting from the point process model, proposed and analyzed in Refs. [9|, I44l - t48j . the nonlinear SDEs generating 
processes with 1//^ noise are derived 

^M,M- The general expression for the proposed class of Ito SDEs is 



dx 



V 



A I x^'^-^dt + ax'^dW . 



(5) 



Here x is the signal, 77 1 is the power-law exponent of the multiplicative noise, A defines the behavior of stationary 
probability distribution, and is a standard Wiener process (the Brownian motion) . 

The nonlinear SDE ([5]) has the simplest form of the multiplicative noise term, ax'^dW. Equations with multiplicative 
noise and with the drift coefhcient proportional to the Stratonovich drift correction for transformation from the 
Stratonovich to the Ito stochastic equation , generate signals with the power-law distributions [l0| . Equation (O 
is of such a type. Therefore, the relationship between the exponents in the drift term, 277 — 1, and in the noise term, 
77, of these equations follows from the requirement of modeling the signals with the power-law distributions. More 
reasoningof the correlation of these exponents, and of the type of equations like ([5]) in general, have been given in 

Refs. I, Hm, [13, mill. 

On the other hand, the simple transformation of the variable, y = x", gives equation of the same type ([5]) only 
with different parameters, a' — aa, rj' = (ji — l)/a -|- 1, and A' = (A — I)/ a + 1. E.g., for a = 1 — 77 with 77 7^ 1 we 
get 77' = 0, i.e., equation for the variable y having additive noise and nonlinear drift, the well-known in econophysics 
and finance SDE describing the Bessel process [43|. Thus the observable x may be a function of another variable y, 
described by a simpler SDE with additive noise. 

Nonlinear SDE, corresponding to a particular case of ([5|) with 77 = 0, i.e., with linear noise and non-linear drift, was 
considered in Ref. (s^ . It has been found that if the damping decrease with increasing \x\, then the solution of such 
a nonlinear SDE has long correlation time. The connection of the power spectral density of the signal generated by 
SDE ([SJ with the behavior of the eigenvalues of the corresponding Fokker-Planck equation was analyzed in Ref. [55| . 
This connection was generalized in Ref. |56| where it has been shown that 1//'^ noise is equivalent to a Markovian 
eigenstructure power relation. 

In order to obtain a stationary process and avoid the divergence of steady state PDF the diffusion of stochastic 
variable x should be restricted at least from the side of small values, or equation ([5|) should be modified. The Fokker- 
Planck equation corresponding to SDE ([5]) with restrictions of diffusion of stochastic variable x gives the power-law 
steady state PDF 



Pix) 



(6) 



with the exponent A, when the variable x is far from the ends of the diffusion interval. The simplest choice of the 



restriction is the refiective boundary conditions at x 
the steady state PDF 



cCjriin and X 



Pix) 



1 



■ exp 



Exponentially restricted diffusion with 



(7) 



is generated by the SDE 



dx — 



1 , 77T, 

77 Ah 

' 2 2 



v^^-'dt 



(8) 
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FIG. 1. (Color online) (a) Steady state PDF P{x) of the signal generated by Eq. (fTO|) with Xmin ~ 1 and Xmux ~ 1000. The 
dashed (green) line is the analytical expression for the steady state PDF. (b) Power spectral density S{f) of the same signal. 
The dashed (green) line shows the slope 1//. 



obtained from Eq. ([5]) by introducing the additional terms. 

In Refs. [13] it was shown that SDE ([5]) generates signals with power spectral density 



S{f) 



1 

7^' 



/3 



1 



A~3 
2(^7-1) 



(9) 



in a wide interval of frequencies. SDE ^ exhibits the following scaling property: changing the stochastic variable 
from a; to a scaled variable x' — ax changes the time-scale of the equation to t' = a^^^^^H leaving the form of the 
equation unchanged. This scaling property is one of the reasons for the appearance of the 1/ power spectral density. 

For A = 3 we get that P — 1 and SDE (O should give signal exhibiting 1// noise. One example of the equation (jH]) 
with A = 3, TO = 1, (T = 1, and rj = 5/2 is 



dx 



1 



1 



x'^dt + x^l'^dW . 



(10) 



Note, that rj = b/2 corresponds to the simplest point process model with the Brownian motion of the interevent time 
Tfe in the events space (fc-space) [3, [H, [s^l , 

dTk = adWk . (11) 

Consequently, the simple point process model may provide one possible reasoning of use of the strongly nonlinear 
multiplicative SDEs for modeling of long-range correlated systems. 

Comparison of numerically obtained steady state PDF and power spectral density with analytical expressions is 
presented in Fig.[T] For the numerical solution we use Euler-Marujama approximation with variable step of integration, 
transforming the differential equations to the difference equations [l^, [50[ . We see a good agreement of the numerical 
results with the analytical expressions. Numerical solution of the equations confirms the presence of the frequency 
region for which the power spectral density has dependence. The width of this region can be increased by 

increasing the ratio between the minimum and the maximum values of the stochastic variable x. In addition, the 
region in the power spectral density with the power-law behavior depends on the exponent rj: if 77 = 1 then this width 
is zero; the width increases with increasing the difference |?7 — 1| |55| . 

The numerical analysis of the proposed SDE ([5]) reveals the secondary structure of the signal composed of peaks or 
bursts, corresponding to the large deviations of the variable x from the proper average fluctuations [13]. Bursts are 
characterized by power-law distributions of burst size, burst duration, and interburst time. 



III. STOCHASTIC DIFFERENTIAL EQUATIONS GIVING g-DISTRIBUTIONS 

The power spectral density of the form l/f^ is determined mainly by power- law behavior of the coefficients of 
SDEs (O, dS]) at large values of x S> Xmin- Changing the coefficients at small x, the spectrum preserves the power-law 
behavior. In addition, Fokker-Planck equation corresponding to SDE ([5]) gives the steady state PDF with power- 
law dependence on x as does the q-exponential function for large x. Therefore, SDE ((S]) can be modified to yield 
generalized canonical distributions of nonextensive statistical mechanics. 
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FIG. 2. (Color online) (a) Steady state PDF P{x) of the signal generated by Eq. (|12p . The dashed (green) line is the analytical 
g-exponential expression (|13p for the steady state PDF. (b) Power spectral density S{f) of the same signal. The dashed (green) 
line shows the slope 1//. The parameters used are A = 3, = 5/2, xo = 1, and cr = 1. 



A. g-exponential distribution 



Modified stochastic differential equation 

dx = a^(^r]^ ^A^ (x + Xof'^'^dt + a{x + Xo)'^dW (12) 

with the reflective boundary condition at cc = was considered in [l^] ■ The Fokker-Planck equation corresponding to 
SDE (IT^ for X > gives g-exponential steady state PDF 

P{x)^^(^^] =^exp^{-Xx/xo), q = l + l/X. (13) 
Xo \x + xqJ Xo ^ 

The addition of parameter Xo restricts the divergence of the power-law distribution of a; at a; —> 0. Equation (1121) for 
small X <ti Xo represents the linear additive stochastic process generating the Brownian motion with the steady drift, 
while for x ^ xq it reduces to the multiplicative SDE ([S]). This modification of the SDE retains the frequency region 
with 1//^ behavior of the power spectral density. 

Comparison of numerically obtained steady state PDF and power spectral density with analytical expressions is 
presented in Fig. [21 We see a good agreement of the numerical results with the analytical expressions. Numerical 
solution confirms the presence of the frequency region where the power spectral density has 1/ dependence. The 
lower bound of this frequency region depends on the parameter xo- 



B. g-Gaussian distribution 

Stochastic differential equation 

dx = a^(^v- {x^ + xlf-^dt + g{x^ + xlY^'^dW (14) 

in contrast to all other equations analyzed in this article, allows negative values of x. This equation was introduced in 
Refs. [57l - [59l |. The simple case 77 = 1 is used in the model of return in Ref. [s^l. Note, that 77 = 1 does not give 1//^ 
power spectral density. The Fokker-Planck equation corresponding to SDE gives q-Gaussian steady state PDF 

g = 1 + 2/A. 

The addition of parameter xo restricts the divergence of the power- law distribution of a; at a; 0. Equation (fT4|) 
for small \x\ <C xo represents the linear additive stochastic process generating the Brownian motion with the linear 
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P(x) 




FIG. 3. (Color online) (a) Steady state PDF P{x) of the signal generated by Eq. (|14|l . The dashed (green) line is the analytical 
g-Gaussian expression (I15|l for the steady state PDF. (b) Power spectral density S{f) of the same signal. The dashed (green) 
line shows the slope 1//. The parameters used are A = 3, = 5/2, xq = 1, and cr = 1. 

relaxation, while for x ^ xq it reduces to the multiplicative SDE ([5]). This modification of the SDE, even the 
introduction of negative values of the stochastic variable x, does not destroy the frequency region with 1/ behavior 
of the power spectral density. 

Comparison of numerically obtained steady state PDF and power spectral density with analytical expressions is 
presented in Fig. |3l A good agreement of the numerical results with the analytical expressions is found. Numerical 
solution confirms the presence of the frequency region where the power spectral density has 1//^ dependence. 



IV. SUPERSTATISTICS AND l/f^ NOISE 



Many non-equilibrium systems exhibit spatial or temporal fluctuations of some parameter. There are two time 
scales: the scale in which the dynamics is able to reach a stationary state, and the scale at which the fluctuating 
parameter evolves. A particular case is when the time needed for the system to reach stationarity is much smaller than 
the scale at which the fluctuating parameter changes. In the long-term, the non-equilibrium system is described by the 
superposition of different local dynamics at different time intervals, that has been called superstatistics ^24. 60-63.]. 
The superstatistical framework has successfully been app lied on a widespread of problems like: interactions between 
hadrons from cosmic rays f25| , fluid turbulence [l^, |64| - [66| , granular material [i^l , electronics [1^ , and economics 

In this article we will consider the case when the fluctuating parameter x evolves according to earlier introduced 
SDE ([5]). The parameter x changes slowly and can by taken as a constant through a period of time T. Due to 
the scaling properties of Eq. ([5]), mentioned in Sec. [TTl the characteristic time scale in Eq. ([5]) decreases as a power 
of x. In order to avoid short time scales and rapid changes of the parameter x, the possible values of x should be 
restricted from above. If the maximum value of the parameter x is Xmax, then the time T during which the parameter 
X changes slowly decreases with increase of Xmax- Within time scale T the signal x has local stationary PDF ip{x\x). 
The long-term stationary PDF of the signal x is determined as 



P{x) 



ip{x\x)p{x)dx . 



(16) 



We can expect that at small frequencies lu <C the spectrum of the signal x is determined mainly by the driving 
SDE. Therefore, we can get the distribution P{x) determined by Eq. (|16p and 1//'^ power spectral density in a wide 
region of frequencies. Using the superstatistical approach, from SDE ([5]) with the exponential restriction of diffusion, 
it is possible to obtain the Tsallis probability distributions. 



A. g-exponential distribution 



In order to obtain g-exponential long-term PDF of the signal x we will consider the local stationary PDF conditioned 
to value of the parameter x in the form of exponential distribution 



ip{x\x) = X ^exp{—x/x). 



(17) 
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FIG. 4. (Color online) (a) Long-term PDF P{x) of the signal generated by Eqs. (jlZp and (jlSp . The dashed (green) line is the 
analytical expression ((20]) for the long-term PDF. (b) Power spectral density S{f) of the same signal. The dashed (green) line 
shows the slope 1//. The parameters used are A = 3, ?7 = 5/2, xo = 1, o- — 1, and Xmax = 10^. 



A Poissonian-like process with slowly diffusing time-dependent average interevent time was considered in Ref. [5S 
The mean x of the distribution ip{x\x) obeys SDE with exponential restriction of diffusion, 



dx 



A 1 xo 1 X 



2 2 X 2 Xniax 



v'^'^-^dt + ax'^dW . 



(18) 



Here a;o is a parameter describing exponential cut-off of the steady state PDF of x at small values of x and the 
parameter Xmax 3> xq leads to exponential cut-off at large values of x. When x ^ Xmax the influence of the 
exponential cut-off at large values of x is small. Neglecting Xmax the steady state PDF from the Fokker-Planck 
equation corresponding to Eq. (jl8|) is 



p{x) 



1 /xo\^ f 



xoT{X - 1) 



(19) 



Using Eqs. (fT6|). (fT7|) . and (fT9|) . we get that for x <C Xmax 
function, 



the long-term stationary PDF of signal x is g-exponential 



P{x) 



A~ 1 

Xo 



Xo 



X + Xo 



A- 1 
Xo 



exp (-Ax/xo) , (7 = 1 + 1/A. 



(20) 



Comparison of numerically obtained long-term PDF and power spectral density with analytical expressions is 
presented in Fig. 2] Numerical solution confirms the presence of the frequency region where the power spectral 
density has dependence. In addition, the long-term PDF of the signal deviates from the g-exponential function 
pi)) only slightly. 



B. g-Gaussian distribution 



In order to obtain the g-Gaussian long-term PDF of the signal x we will consider the local stationary PDF condi- 
tioned to the value of the parameter x in form of the Gaussian distribution. 



ip{x\x) 



1 



■ exp(— x^/a;^) . 



(21) 



The standard deviation of x in the distribution ip{x\x) is proportional to the parameter x. The fluctuating parameter 
X obeys SDE with exponential restriction of diffusion ([8]) with the parameter to = 2, 



dx 



v- 



X 
2 



i^'^-^dt + ax'^dW . 



(22) 



Here xq is the parameter describing exponential cut-off of the steady state PDF of x at small values of x, whereas 
the parameter a;max(^ xq) leads to the exponential cut-off at large values of x. When x <C imax, the influence of 
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FIG. 5. (Color online) (a) Long-term PDF P{x) of the signal generated by Eqs. pip and H22p . The dashed (green) line is the 
analytical expression 1)24^ for the long-term PDF. (b) Power spectral density of the same signal. The dashed (green) line shows 
the slope 1//. The parameters used are A = 3, = 5/2, xo = 1, a — I, and imax ~ 10^. 



the exponential cut-off at large values of x is small. Neglecting Xmax the steady state PDF from the Fokker-Planck 
equation corresponding to Eq. is 

From Eqs. ([T6|. ((2T|) . and pB]) we obtain that for x ^ Xmax the long-term stationary PDF of the signal x is q-Gaussian, 
i.e., 



Pix) ^(^) " - ^(^) exp (-X^) (24) 

q=l + 2/X. 

Comparison of numerically obtained long-term PDF and power spectral density with analytical expressions is 
presented in Fig. [5] Numerical solution confirms the presence of the frequency region where the power spectral 
density has 1//^ dependence. In addition, the long-term PDF of the signal deviates only slightly from the g-Gaussian 
function In contrast to Sec. IIIIBl the superstatistical approach yields 1//^ power spectral density only for the 

absolute value of the signal. Since the signs of two consecutive values of x are uncorrelated, the spectrum of the 
signal X itself in the same frequency region is almost flat. 



V. DISCUSSION 



Common characteristics of complex systems include long-range interactions, long-range correlations, (multi)fractality, 
and non-Gaussian distributions with asymptotic power-law behavior. The nonextensive statistical mechanics, which 
is a generalization of the Boltzmann-Gibbs theory, is a possible theoretical framework for describing these systems. 
However, long-range correlations and noise has not been observed in previously used models giving distributions 
of the nonextensive statistical mechanics. The joint reproduction of the distributions of nonextensive statistical 
mechanics and 1/f noise, presented in this paper, extends understanding of the complex systems. 

Modeling of the concrete systems by the nonlinear SDEs is not the goal of this paper. However, relation of the models 
and obtained results with some features of the financial systems can be pointed out. Equations with multiplicative 
noise are already used for modeling of financial systems, e.g. the well-known 3/2 model of stochastic volatility (75| . 
There is empirical evidence that trading activity, trading volume, and volatility are stochastic variables with the 
long-range correlation i42, 43, 75-78]. This key aspect, however, is not accounted for in the widely used models. 
On the other hand, the empirical findings of PDF of the return and other financial variables are successfully described 
within the nonextensive statistical framework p^ . [28l [29j . The return has a distribution that is very well fitted by 
g-Gaussians, only slowly becoming Gaussian as the time scale approaches months, years and longer time horizons. 
Another interesting statistic which can be modeled within the nonextensive framework is the distribution of volumes, 
defined as the number of shares traded. Modeling of some properties of the financial systems using the point process 
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models and SDEs have been undertaken in papers [52, l58|, [591, l79| . Equations presented in this article incorporating 
the long-range correlations, 1/ noise, and g-Gaussian distributions suggest deeper comprehension of these processes. 
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